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We present a new technique that improves the scaUng of the error in the adiabatic approximation 
with respect to the evolution duration, thereby permitting faster transfer at a fixed error tolerance. 
Our method is conceptually different from previously proposed techniques: it exploits a commonly 
overlooked phase interference effect that occurs predictably at specific evolution times, suppressing 
transitions away from the adiabatically transferred eigenstate. Our method can be used in concert 
with existing adiabatic optimization techniques, such as local adiabatic evolutions or boundary can- 
cellation methods. We perform a full error analysis of our phase interference method along with 
existing boundary cancellation techniques and show a tradeoff between error-scaling and experimen- 
tal precision. We illustrate these findings using two examples, showing improved error-scaling for 
an adiabatic search algorithm and a tunable two-qubit quantum logic gate. 
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The adiabatic approximation underpins many important present-day and future applications, such as stimulated 
rapid adiabatic passage (STIRAP) coherent control of chemical reactions Q, and quantum information process- 

ing (QIP) This approximation asserts that a system will remain in an instantaneous eigenstate of a time- varying 

Hamiltonian if the time- variation happens slowly enough. Errors in this approximation correspond to transitions away 
from the instantaneous ("adiabatically transferred") eigenstate. For high-performance applications, it is not always 
practical to minimize errors by slowing things down. Ambitious future technologies, such as quantum computing 
—I . devices, will demand simultaneous maximization of both accuracy and speed. 

I In this paper, we investigate a phase cancellation effect that appears during an adiabatic evolution and can be 
^ , exploited to polynomially reduce the probability of a given transition at fixed maximum evolution time. This can 
O^' lead to speed increases at fixed error probability. Unlike alternative methods that obtain improvements by modifying 
the adiabatic path 0], our technique chooses the evolution time so that destructive interference suppresses the 
' transition. Furthermore, this phase cancellation effect can be exploited to improve existing adiabatic error reduction 
I strategies such as local adiabatic evolutions or boundary cancellation methods. We provide an error analysis of our 
00 ■ method and conclude that the accuracy improvements come at the price of increasingly precise knowledge of the time- 
' dependent Hamiltonian; this implies that accuracy is an important and quantifiable resource for quantum protocols 
utilizing adiabatic passage. 

in 

O ■ I. ADIABATIC APPROXIMATION 

Following previous authors [1, Q, we define the error £ to be the component of the post-evolution state vector 
that is orthogonal to the state intended for adiabatic transfer. For convenience, we represent all mathematical terms 
as explicit functions of the "reduced time" s{t) = t/T, where t is the time, T is the total evolution duration, and 
! < s < 1. This parameterization leaves the form of the Hamiltonian 7^(s) unchanged as T varies. 
■ In many circumstances, the following criterion adequately estimates the magnitude of the total error £ at time 
t = T for a given Hamiltonian: 

,|.,I<1 11^^(^)11 

f <-max^ r^T~\ f ( M2 ' 1) 

T s mm^ I (s) - So (s)r 

where Ei,{s) {v ^ 0) is the instantaneous energy of the v^^^ eigenstate of the Hamiltonian 'H(s) and Eq{s) is the energy 
of the eigenstate being transferred (usually the ground state) 0, We use the convention h= 1. 

Although cq. ^ provides an expedient heuristic for estimating the accuracy of adiabatic passage, it is (in general) 
neither necessary nor sufficient to bound the fidelity of adiabatic state transfer [ll|, [T^] ■ This equivocality opens the 
possibility of a modest allocation of resources being used to enable significantly improved error-scaling. 

One method of improving the fidelity of adiabatic transfer is via the use of a "local adiabatic" evolution d, 0, [ll] ■ 
The idea behind the local adiabatic approximation is to tailor variation of H with respect to s to minimize the 
instantaneous nonadiabatic transition rate \\-^'H{s)\\/i[mn,^\E^{s) — Eq{s)\'^ . Local adiabatic methods have lead to 
substantial improvements in the asymptotic error-scaling 8 with respect to the Hilbcrt space dimension N 0, [l3| ; 
however, they do not improve the scaling of the error with T. 
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Recently, methods were developed for improving the scaling of £ with T from order 0{1/T) to 0(1/T™+^) by 
setting the first m derivatives of the Hamiltonian to zero at the beginning and end of the evolution [14 . [l5| . Error 
reduction techniques employing this result are collectively referred to as "boundary cancellation methods." Boundary 
cancellation methods have two main drawbacks: first, they assume that the first m derivatives of 'H(s) are exactly 
zero, leaving it unclear whether they are robust against small variations in the derivatives of the Hamiltonian; second, 
in the regime of short T these methods can have error-scaling that is inferior to the trivial case wherein no boundary 
cancellation technique is applied (i.e., to = 0). Our work addresses these problems: we first provide an analysis of 
the sensitivity of boundary cancellation methods to small variations in the values of the first to derivatives of 'H(s); 
we then show that phase interference can be used to further reduce errors without increasing m, providing better 
error-scaling for short T. 



II. MAIN RESULT 



We present a new technique for quadratically suppressing the probability of a particular nonadiabatic transi- 
tion during adiabatic passage. It works by exploiting a phase interference effect that appears in adiabatic systems 
with Hamiltonians obeying a simple symmetry. This effect can be exploited in a realistic class of time-dependent 
Hamiltonians that includes many adiabatic algorithms and transport protocols, as well as any Hamiltonian obeying 
H{Q)=H{1). 

Consider a time-dependent Hamiltonian 'H(s) acting on an ^''-dimensional Hilbert space spanned by the instanta- 
neous energy eigenvectors \v{s)) where = 0, 1, . . . , iV — 1. We define |0(s)) to be the state intended for adiabatic 
passage. We use the notation T-L''P^{s) = {-§^yH{x)\s. In Section Hh] we wiU show that errors in adiabatic passage can 
be reduced for Hamiltonians obeying the boundary symmetry condition. 

Ml)|H('"+i)(l)|0(l)) _ ( (KO)|H("+i)(0)|0(0))^ ^2) 



where 9 is an arbitrary phase factor, and to is the number of derivatives oi'H{s) that are zero at the boundaries s = 0, 1 
(e.g., if TO = 2 then the first and second derivatives of T-L{s) are zero at the boundaries, whereas if 771 = then none 
are zero on the boundary). In practice, any time-dependent Hamiltonian may be adapted to satisfy cq. ([2]), simply 
by adjusting its rate of change in s at the boundaries. If eq. ^ is not exactly satisfied then the phase interference 
effect will still reduce errors, but it will not necessarily improve the asymptotic error-scaling with T. 

Our method can be used in conjunction with existing boundary cancellation methods to produce even greater 
improvements in the asymptotic error-scaling with T. Amplitudes of the transitions |0(0)) — ?> |j^(l)) are reduced from 
the order 0{T~™'~^) estimates given in refs. [13, [3 to order ©(T"™"^) at the discrete set of times T = Tn,^, where 
n is an even integer and 

rnr — 6 

This can lead to polynomial reductions in the overall error if \\£\\ is dominated by a small number of transitions. 

We refer to boundary cancellation methods that are augmented by our scheme to produce order 0{T~"^~'^) error- 
scaling as "augmentcnted boundary cancellation methods." In Section IIVI we will analyze the error robustness of 
our augmented boundary cancellation method along with the original schemes laid out in refs. [13 ■ We show 
that performance improvements are derived from accurate knowledge of the system's eigenspectrum {E^}, its total 
evolution time, and the derivatives of its Hamiltonian, and we provide quantitative error-bounds on these quantities. 
We provide numerical examples that verify the predictions of our theory in Sections |V] and IVII 



III. THEORY 



We will break our discussion of the theory of our method into two parts. First, we discuss the special case for which 
m = 0. This simple case is conceptually distinct from existing boundary cancellation techniques, which require to > 
to produce improvements over eq. (P). We then discuss the more general case in which to > 0. 

To obtain our results, it is not necessary to assume that the instantaneous eigenvalues satisfy the ordering condition 
Eo{s) <Ei{s) < . . .<i?jv-i(s)- We do however require that Eo{s) ^ E^{s) V J/ > 0, unless transitions between |0(s)) 
and |7^(s)) are strictly forbidden by H(s). For convenience, we also assume that the phases of the instantaneous 
eigenvectors are chosen such that (t'(s)|j/(s)) = 0. This choice does not affect the quantum dynamics, but it simplifies 
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the analysis of the error. We also assume that the Hamiltonian is difFerentiable m + 2 times and that each derivative 
is bounded for all T . These last restrictions are put in place order to prevent issues that arise for Hamiltonians 
resembling that of the Marzlin-Sanders counterexample [ill, [iBl ■ 

Given the above assumptions, the error in the adiabatic approximation £ for a Hamiltonian evolution acting on an 
A'^-dimensional Hilbert space is given by 



Af-l 



(4) 



i/=i 



We know from previous work that £u € 0{T~"^~^) [ij, [1^, and asymptotically tight expressions are known for £i, 
in the to = case [1, [lB| • We therefore begin with this case to illustrate how our phase interference effect can be 
utilized. Given that to = 0, the form of £i, reduces to 



£. 



(5) 



and where '-f^{s)~ Eo{s) — Ei,{s). If wc choose H(s) to obey then the absolute value of eq. ^ reduces to 



(6) 



Eq. (j6l) has extrema at T = Tn,u- It is maximized when n is odd and vanishes when n is even. Thus, when T = Tn.u 
(even n), phase interference causes the magnitude of the z^'^ component of £ to be quadratically reduced from ©(T^^) 
to 0(T-2). 

If TO > then the phase interference effect also suppresses probability of excitation to \v{l)) at T = T„.iy for any 
even integer n > 0, but this effect does not directly follow from existing results. We show in Appendix 1X1 using a 
perturbative expansion (similar in reasoning to that of refs. [T^. [T5|) that if the first to derivatives of 'H(s) are zero at 
the boundaries s = 0,1 then 



\£.\ = 



(i^(s)|-H(™+i) (s)|0(s))e-* /o 7.(e)d«T 



s=0 



(7) 



Similar to eq. (O, eq. ([7]) reveals an adiabatic phase interference effect also that suppresses the error at certain 
times. This suppression occurs when 



(i^(l)|H(™+iHl)|0(l))e-'^^oN.(^)d. (z.(0)|H(™+iHO)|0(0)) 



7. (!)"+' 



7.(0) 



m+2 



(8) 



implying that adiabatic phase interference effects reduce the order of transition amphtude from 0{T 



-rn— 1 \ 



to 



0{T 



'm-2 



) when T — Tn^i, for even n. 



IV. TOLERANCES 



Limits on the precision of physical apparatus prevent perfect phase cancellation in realistic applications. Errors 
can result from imperfect modelling of the Hamiltonian, inexact calculations of the gap integrals, or inaccuracies in 
the timing or control apparatus. It is therefore necessary to address the impact of empirical imperfections on the 
feasibility of augmented boundary cancellation methods and determine when they methods can be experimentally 
realized. 

"Symmetry errors" occur when the timing symmetry condition ^ is not precisely satisfied: 



(Kl)|H(™+i)(l)|0(l)) (KO)|H(™+i)(0)|0(0)) , 

— e 



7.(1) 



m+2 



7.(0) 



m+2 



> 0. 



(9) 



Comparing eq. ^ with eq. (O, we find that the contributions to £^ due symmetry errors are of order 0{T' 
long as AS"^ e 0{T-^). 



-m — 2\ 



4 



"Gap errors" occur when inaccuracies in the estimate of the gap integral leave condition ^ unsatisfied: 



AG. 



> 0. (10) 



Expanding eq. (O in powers of AG,y, we find that the contributions to due to gap errors arc of order 0{T ^) 
if AG. e 0{T-^). 

"Timing errors" occur when the actual evolution time T differs from the ideal evolution time T„_.: 

Ar„,, = |T„,, - T| > 0. (11) 

Expanding eq. (O in powers of AT„ we find that the contributions to 5. due to timing errors are of order ©(T""'^^) 
if AT„,. eO(T-i). 

"Derivative errors" can also occur wherein one or more of the derivatives of the Hamiltonian that is assumed to be 
zero is not: 

AH(p) = max 11^(^^5)11 > 0, (12) 

s— 0,1 

for p ~ 1, . . . , m. Such errors do not affect the error-scaling if for all such p, 

An^P^ G 0{l/T"'+^-P). (13) 

In other words, given that the first m derivatives of H are approximately zero at the boundaries, the uncertainty in 
each derivative must shrink polynomially as T increases in order to achieve the full promise of an augmented boundary 
cancellation method. The proof that this criteria is sufficient is not simple: it requires a high-order perturbative 
analysis of the error in the adiabatic approximation. Details are provided in Appendix [Bl 

If m is a constant, then it follows that augmented boundary cancellation methods are error robust in the sense 
that their error tolerances scale polynomially with T^^. This is not problematic for numerical studies because 
additional precision can be provided at poly-logarithmic cost. However, experimental errors cannot always be so 
conveniently reduced, and boundary cancellation techniques that use a large value of m may be impractical. The 
situation is even worse if exponential error-scaling is required, which can be obtained if to e 0(r/logT). In such 
circumstances the tolerances ^(^•'(s) decrease exponentially with T and therefore boundary cancellation methods are 
not error robust. This implies that boundary cancellation techniques (augmented or not) cannot in practice achieve 
exponential scaling without exceedingly precise knowledge of the derivatives of the Hamiltonian at the boundaries. 
The m = method may therefore be more experimentally relevant than its higher-order brethren, because of its 
minimal precision requirements and its superior scaling for modestly short T. 

As the performance improvements provided by boundary cancellation methods come at the price of increasingly ac- 
curate information about the Hamiltonian and the evolution time, such information may be viewed as a computational 
resource for protocols utilizing quantum adiabatic passage. This suggests that current analyses [l^l of the resources 
required for generic adiabatic quantum computing may be incomplete. We illustrate this subtlety in Section |V] by 
showing how to quadratically improve the total error-scaling \ \£\\ of an already "optimal" quantum algorithm. 

V. SEARCH HAMILTONIANS 

Adiabatic quantum computing (AQC) algorithms arc natural candidates for error suppression by our technique. 
To demonstrate, we examine an algorithm that adiabatically transforms an initial guessed state into the sought state 
of a search problem [l^] . The Hamiltonian for this algorithm is 

n{s) = / - (1 - <^(s))|+®")(+®"| - (/.(s)|0®")(0®"|, (14) 

where |+) = (|0) + \l))/\/2, |0®") is the state that the algorithm seeks, and (j) : [0,1] ^ [0,1] obeys (j){0) = and 
<^(1) = 1. 

Two common choices for (j){s) 0, [13; 13 't'{s) — s and 

VA^- 1 - tan [ arctan(\/A^- 1)(1 - 2s)] 
s = ^ ; — ^. 15) 
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FIG. 1. Final error amplitude \£\ as a function of T for 
the Search Hamiltonian (|14[l using A*" = 16 and <j){s) — s. 



FIG. 2. Final error amplitude |£^| as a function of T for 
the Search Hamiltonian (|14p using A'^ = 16 and eq. (|15p . 



The latter choice ()15p is said to generate a "local" adiabatic evolution 0, In each case, the (dimensionless) energy 
gap is 



7i(s) = ^1 - 4 (^1 - 1^ <t>{s){l - (16) 

where |0(s)) is the ground state of eq. ([Ti]) and |l(s)) is the only other eigenstate that is coupled to |0(s)) [l^. From 
the eigenvectors of H{s), it is straightforward to verify that both forms of (/)(s) given above satisfy eq. with m = 0. 

Figs. [1] and [2] show that the choice T = Tn^v (even n) produces quadratic improvements in the scaling of for 
both (/)(s) = s and eq. (jlSp at large T . For odd values of n, the error is maximized, as expected. It is apparent that 
randomly selected times are extremely unlikely to exhibit maximum phase cancellation. Figs. [1] and [2] also suggest 
a second benefit of our technique: existing boundary cancellation methods [3, [13 '^^"^ improve the performance of 
adiabatic algorithms in the limit of large T, but these improvements come at the price of inferior error-scaling for 
small T . as seen in Fig.[3]of ref. [H. The results shown here in Figs. [1] and [2] exhibit no such tradeoff. 

Figs. [Hand [5] also shed light on the nature of the complexity of adiabatic algorithms. Several previous studies have 
taken the complexity of an adiabatic algorithm to be given by the evolution time required for the error predicted by 
eq. ID) to fall within a specified tolerance 0, [13: [13 ■ I^^ '^^^'^ of the local adiabatic evolution, this time scales 
as 0(\/]V), which is known to be optimal 0, [l3l- Fig- El show that this error can still be quadratically reduced by 
eliminating the 0(T~^) contributions to it. These results do not violate quantum lower bounds because the time 
required for the 0{1/T) to become dominant still scales as 0{\/N) [i3|. Therefore even an exponential improvement 
in the subsequent adiabatic regime would not violate quantum lower bounds. Paradoxically, these results suggest that 
the complexity of adiabatic algorithms may be dictated by the physics of the sudden approximation rather than the 
adiabatic approximation. 

We demonstrate our generalized m > technique in Fig. [3l where we plot jf^j as a function of the total evolution 
time for a search Hamiltonian with 0(s) taken to be 

f'a;"(l-a;)™dx , ^ 

/q x™(1 - x)"da; 

This interpolation was originally suggested in ref. [31 and is chosen because it conveniently guarantees that the first 
m derivatives of 'H(s) are zero at s = and s = 1. Additionally, in the m = case it gives the linear interpolation 
(j){s) = s used in Fig. [TJ 

Fig. 13] demonstrates the improvements that arise from combining our results with those taken from refs. [T^. [isj. 
It is notable to see that the m = data in Fig. [3] (a) nearly coincides with that for m = 1 in Fig. [3] (b) for sufficiently 
large T. Similarly, the m = 1 data in Fig. [3] (a) corresponds to the m = 2 data in [3] (b) in the same limit. This 
shows that our technique can be used to improve the overall accuracy of boundary cancellation techniques without 
compromising the error-scaling for short T . 



6 




1000 2000 3000 4000 5000 6000 7000 8000 1000 2000 3000 4000 5000 6000 7000 8000 

T T 

(a) Cb) 

FIG. 3. This figure shows that existing boundary cancellation methods can be augmented with our boundary cancellation 
method to achieve even higher-order error-scaling for a search Hamiltonian with A'^ = 16 and 0(s) chosen as in eq. (|17|l . Fig. |3] 
(a) is a plot of the error at the times when our theory predicts improved error-scaling (i.e., even n), whereas Fig. [3](b) displays 
the times when the errors are predicted to be maximized (i.e., odd n). 



VI. TWO-QUBIT GATE. 



Our technique naturally lends itself to Haniiltonians that couple the ground state to only one excited state, such 
as the Search Hamiltonian given in eq. (jl4p . If the total error \\£\\ is dominated by several transitions, this technique 
can still be adapted to approximately cancel multiple transitions simultaneously. To demonstrate, we show how to 
optimize the fidelity of an adiabatic two-qubit logic gate without decreasing its speed. Similar improvements were 
reported previously poj . without a broadly-applicable underlying theory or error bounds. 

We apply of our method to an exchange-based two-qubit operation designed for neutral atom QIP 0, [2ll - [23l | . This 
operation exploits identical particle exchange to generate a partial "swap" operation between qubits stored in nuclear 
spin [22I or valence electronic states 2^ of optically trapped atoms. The gate generates a relative phase of e~*" 
between the symmetric and antisymmetric components of the particles' vibrational degrees of freedom. The phase 
difference is then transferred to the respective components of the two-qubit subspace '■ i,j G {0,1}}. This 

produces an operation that (with single-qubit rotations) is locally equivalent to a tunable entangling controUed-phase 
gate e-2-|ii><ii|[2|. 

Following previous work [H, [2^ , we examine a simple Hamiltonian governing two identical particles confined to one 
dimension and trapped by pair of moving potential wells. The Hamiltonian for particles 1 and 2 is given by 



'H{xi,X2,Pi,P2,s) = 'H{xi,pi,s) + 'H{x2,P2,s) + 2auj±S{xi-'X2), (18) 

for ^{{XjP, s) = p^ I2m + V{x+ {s — \ )d) + y (x — (s— \ )d), where x and p are the position and momentum of a particle 
of mass m. The potential V{x) ~ —Vo exp(— x^/2cr^) describes a ID Gaussian trap of depth Vo and variance . Traps 
are initially separated by a distance d = Sc. We consider a ID s-wave scattering interaction, with scattering length 
auj± = 3(7 and transverse confinement frequency uj± [24| . As eq. (jl8p is symmetric, transitions between symmetric 
and antisymmetric states are forbidden, and each symmetry subspace evolves independently. 

We diagonalized eq. ([T8|) over the range < s < 0.5 at As = 1/1200 intervals. We then used a spline fitting to 
integrate eq. obtaining numerical estimates T of the ideal r„ The quality of initial approximations were then 
improved using the relationship |T„^,y — Tj « T/An, where An measures the beat frequency between and T^l 
(e.g., the distance between cusps on inset, Fig.|4]). More sophisticated model Hamiltonians may be solved using more 
advanced numerical techniques and empirically refined in the same manner. 

We numerically integrated the Shrodinger equation to obtain system dynamics of durations {T^ 5}, explicitly gen- 
erating sets of wave functions {IV'n (s))} and {|'/'^(s))} for two distinct initial states: the symmetric ground state 
IV'Tt(O)) = |0(0)) and the antisymmetric (effective) ground state IV'tTCO)) = We chose T„_5 because |5(s)) is the 

first eigenstate that significantly couples to |0(s)). This transition is dominant because the Of^l, 0o2, and f-)- 3 
transitions are forbidden, and the 0<-^4 coupling is weak. We define KV^I*^)! = KV'n 
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FIG. 4. Transition amplitudes and bounds for T — T„^5, 
over 200 < n < 2000. Main figure shows \{^pt\5)\ « l^s 
for even n (soUd) and odd n (dashed), which are bounded by 
ma.Xs[2\\-j^H{x,p,s)\\ /{E5{s)-Eo{s))'^] (dotted). Inset shows 

|{V';;|6)| bounded by max,[2 \\^^H{x,p, s)\\ /{E7{s)-Ei{s)f]. predicts \\£+f < 0.046 and \\£-f < 0.62 x 10'^ at n = 460 



TABLE I. Error probabilities and the phase gap a (radians) 
obtained from simulation runs {Tn.s} for 456 < n < 474. For 
these times, local minima of |(i/'n|5)| roughly match those of 
|(i/'^|6)| (inset. Fig. 2]) and |(?/j^|7)|. We denote total errors 
as \\£+f = 1 - \{^tmf and \\£-f = 1 - ■ Eq. © 



The error probabilities are improved by nearly three orders of magnitude over the bound set by cq. ([1} by applying 
our technique to this system (Table |T| . This corresponds to a tenfold increase in gate speed (given a maximum error 
rate of 10~*), for the linear motion described by eq. ([T5)) . Greater improvements could be achieved by choosing 
'H{x,p,s) or s{t) to satisfy eq. ([3]) for more transitions simultaneously and with better synchronization. 

Partial svifap operations have been experimentally demonstrated using neutral atoms in a double- well optical lattice, 
but the adiabatic requirement limits gate times (~4ms for high fidelity operation [2l|). Our technique thus affords a 
significant advancement to inherently slow gates of this kind. Furthermore, because the phase a scales with T (see 
Table |T| , the precision necessary for accurate gate operation is itself comparable to that needed to implement our 
phase cancellation technique on an atomic quantum logic gate. 

We have numerically demonstrated that error in the adiabatic approximation can be reduced for an experimentally 
relevant model of a quantum gate. An important remaining issue is whether the experimental uncertainties required 
to observe error reductions are reasonable for this model system. By first-order Taylor expansion of eq. ([5]) , we find 
that if 

AS, I i < 33% and = #^ < 0.02%, (19) 

/ V 75 (0)7 /o 75(s)ds 7460,5 

then the observed transition amplitude at T sa ?46o,5 will be less than half of that at T = 2459.5. These modest 
requirements imply that our m = method may be rapidly incorporated into present-day or near-future atom-based 
QIP experiments. Such an experiment would also provide a highly sensitive test of the validity of the adiabatic 
approximation in open quantum systems. 



VII. CONCLUSION 



We have presented a new technique for improving the fidelity of adiabatic transport. Our technique exploits 
an adiabatic phase cancellation effect that occurs at certain evolution times to produce improved error-scaling. In 
addition, our method applies directly to a host of experimentally relevant physical systems, often without modification 
to the adiabatic path s{t). Our technique can also be used to improve the accuracy of existing boundary cancellation 
techniques, providing improved scaling over those methods when an easily-satisfiable symmetry condition ([2]) is met. 
We show that these "augmented" boundary cancellation techniques can provide unsurpassed accuracy, requiring 
comparably precise control over the Hamiltonian to achieve high-order error-scaling. Consequently, our work reveals 
that precision (in addition to energy and time) is a subtle and important resource to consider when devising algorithms 
and experiments that utilize adiabatic state transfer. 

We have illustrated these claims using numerical examples of QIP applications. We numerically demonstrated the 
use of augmented boundary cancellation methods for to = 0, 1, 2 for an adiabatic search algorithm. We also optimized 
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a simple adiabatic quantum logic gate using our m = method. In that case we also performed an error analysis and 
found that the error tolerances needed to apply the method are experimentally reasonable. 

Our results open several interesting avenues of further inquiry. We have shown that our technique can be used to 
improve the accuracy of some local adiabatic evolutions, but it would be interesting to see if further improvements 
can be obtained by using our method in concert with more sophisticated adiabatic optimization methods such as the 
one given in ref. |6| . In addition, determining the error tolerances for small deviations along the adiabatic path would 
be an important step towards fully characterizing precision as a resource for adiabatic processes. Our preliminary 
estimates suggest that it may be possible to observe error reductions for atom-based quantum logic using optical 
dipole traps, but other experimental setups may also be well-suited to study this effect, such as nuclear magnetic 
resonance (NMR) systems. Such experiments would not only be interesting as a test of the viability of augmented 
boundary cancellation methods as an error-reduction strategy, but would also provide a highly sensitive test of the 
limits of the adiabatic approximation itself. 



Appendix A: Proof of Eq. \7\ 

In Section imi we claimed that phase cancellation can be used to accelerate the convergence of boundary cancellation 
techniques. Specifically, we claimed that our augmented boundary cancellation methods reduce \£^\ from order 
©(T-™-!) to 0(r-™-2)^ -y^g ^iii j^Q^ justify why this is the case. 

Using the path-integral representation of the time-evolution operator presented in 0, [l^, [l^ we have that 



1^.(1)1 



•Jo „ "'0 JO 



(Al) 



where (3^^^ is defined for any i> and fi in the set {0, . . . , TV — 1} by 



if E,{s) ^ E^{s), 

E^{s)-E^{s) otntiwise. 



(A2) 



We analyze the series under the assumption that the first m derivatives of the Hamiltonian are zero at the boundaries 
s = 0, 1. Using integration by parts, we find that 



Then, using the fact that 'H^^'(O) = '^'■^•'(1) = 0, the first term on the right side of eq. (|A3[) is zero. Evaluating the 
second term using integration by parts, we obtain 



"'O 



1 'd_ 

ds 7i 



-TV^ (1^4%)) e-^/o -(«)<^«-d.. (A4) 



As before, the first term in this expression is zero because the first two derivatives of the Hamiltonian are zero. We 
then continue by this reasoning, applying integration by parts m + 1 times. Then after dropping the first m derivatives 
of the states, Hamiltonian, and energy gaps at s = 0, 1, we find 



1 



7.(0) 



m+2 



+ C'(l/T™+2) 



(A5) 



We then see that the symmetry condition in eq. ^ implies that if T = Tn^u then the first term in the expansion 
in eq. (jASp is 0(l/r™+^). The result of eq. ^ then holds if the remaining terms in eq. (|A3|) are asymptotically 
negligible. 

Turning our attention the remaining path- integrals in eq. (|A3[) . we find that all of the remaining terms are 
0(1/T™^'^). This is because these terms involve contain multiple products of Therefore, if we perform in- 

tegration by parts m -J- 1 times on the outermost integral, then the term involving becomes multiplied by at 
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least one P^^n term, which is zero on the boundary by definition. Therefore, no non-zero terms appear in the expansion 
of these integrals to ©(l/T™"*"^). Hence, the first term in eq. (|A3p is asymptotically dominant as anticipated [sl. [l6j. 
Since the first term is asymptotically dominant and also of order ©(l/T'""*"^) given T = Tn.u-, the augmented boundary 
cancellation technique proposed in Sec. IIIII combines with existing methods. 



Appendix B: Error-Robustness of Augmented Boundary Cancellation Methods 

In Section ITVl we claimed without proof that if the uncertainty in the p^^ derivative of 'H(s) is C)(r^™^^+'') for all 
p = 1, . . . , m, then that derivative can safely be assumed to be negligible. We prove this now by demonstrating that 
the leading order terms involving H^p^(O) or TiSP^l) for p = 1, . . . ,m are of order 0{T~"^~^) under this assumption. 

We begin by assuming that, for some q, H^'^^s) is non-zero at the boundaries s = 0, 1 and that all lower derivatives 
are negligible there. Following the argument put forward in Appendix \^ the lowest order term that appears after 
applying integration by parts q times to eq. (|A1[) is 

1 / (Kl)|H(^)(l)|0(l))e"'/o ^-(-)'^-^ MO)|H(^)(0)|0(0)) \ 

If and n^i^O) are both of order 0(T-'"-2+9), then the term ^ is reduced to order C'(T-'"-2). As 

argued in Appendix \^ other terms that appear in the perturbative series after repeated integrations by parts are 
asymptotically smaller than this term and therefore do not affect the error-scaling. Thus, it is sufficient to render 
errors in the g*^ derivative of 'H(s) negligible by taking them to be 0{T~™~'^^'^). 

By the same reasoning, if the uncertainty in the p*'' derivative of 'H(s) is 0(7"^"^^+^) for all p = 1, . . . , m, then 
the total contribution of derivative errors is 0(T^™^^) given that m £ 0{l). This implies that augmented boundary 
cancellation methods are robust to derivative errors given that to is a fixed integer. This result also trivially implies 
that existing boundary cancellation methods are robust to derivative errors under the same circumstances. 

If TO is not bounded from above by a constant then this analysis fails because the previous analysis ignored 
multiplicative factors of m that appear in the analysis. Such terms could make the neglected higher-order derivative 
terms much larger if to is an increasing function of T. This means that if we wish to achieve exponential error-scaling 
by taking m € 0(r/log(T)), then the tolerance for derivative errors must shrink even further from the already 
exponentially small error tolerances obtained by substituting to € 0(T/log(T)) into 0(T~"^~'^^p) for fixed p. We 
conclude that boundary cancellation methods that exhibit exponential error-scaling are not robust to derivative errors. 
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